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SUMMARY 


A  review  is  given  of  recent  work  on  boundary  layers  in  electrically- 
conducting  fluids.  This  includes  a  discussion  of  the  differences  between 
ordinary  boundary  layers  and  plasma  boundary  layers  with  sp<5cial  reference 
to  the  new  dimensionless  parameters  for  the  magnetic  case,  including  those 
which  arise  from  the  microscopic  structure  of  the  fluid.  The  convection  of 
current  and  vorticity  along  magnetic  field  lines  by  the  Alfven  wave  mechanism 
is  shown  to  be  the  principal  new  propagation  mechanism  in  magnetohydro¬ 
dynamic  boundary  layers.  Recent  work  in  inertial  boundary  layers,  channel 
flows,  wakes,  and  electrode  boundary  layers  is  summarized. 
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I.  Introduction 


While  the  physical  conditions  under  which  a  plasma  may  be  created 
or  maintained  cover  many  orders  of  magnitude  change  in  the  particle  density 
and  temperature,  the  plasma  boundary  layers  which  we  shall  discuss  can 
exist  only  under  the  limited  conditions  for  which  the  plasma  may  be  treated 
as  a  continuum  in  the  usual  fluid'-mechanical  sense.  This  limitation  reflects 
in  part  the  tendency  of  acrodynamicisls  to  move  from  the  familiar  to  the  less 
familiar  by  the  extrapolation  of  existing  understanding  and  in  part  the  possi¬ 
bilities  of  the  practical  application  of  magnetohydrodynamics  (MHD).  The 
plasma  jet,  MHD  generator,  MHD  propulsion  motor,  MHD  flow  meter  and 
other  devices  i\re  already  sufficiently  practical  to  require  a  more  complete 
understanding  of  their  operation,  an  understanding  which  can  probably  be 
acquired  by  the  extrapolation  c{  known  techniques  (such  as  boundary  layer 
theory).  This  should  not  blind  us  to  the  fact  that  only  a  limited  region  of  the 
plasma  state  is  being  subject  to  scientific  and  engineering  scrutiny--a  state 
of  affairs  which  cannot  be  considered  to  be  satisfactory. 

In  the  classical  sense,  a  boundary  layer  is  a  region  of  the  flow  into 
which  disturbances  diffuse  frorp  boundaries,  carrying  mass,  momentum  and 
energy  (including  magnetic  energy).  From  a  practical  point  of  view,  a  boundary 
layer  is  the  region  which  "insulates”  the  flow  from  its  boundaries,  and  in  many 
plasma  devices  the  boundary  layer  provides  the  mechanism  for  containment 
of  the  high  temperature  gases.  In  this  respect,  a  magnetic  field  may  be 
considered  tc  be  a  bounding  surface,  and  therefore,  a  source  of  disturbances 
as  well  as  a  mechanism  for  containment.  It  is  this  point  of  view  which  will  be 
found  most  helpful  in  understanding  MHD  bciuidary  layers. 


In  some  plasma  devices  an  electric  current  must  pass  from  a  solid 
boundary  into  the  fluid.  The  detailed  mechanism  by  which  electrical  charges 
flow  from  the  fluid  to  the  wall  are  not  well  understood  and  appear  to  be  highly 
specific  to  the  wall  properties  and  chemical  state  of  the  plasma.  Thus, 
boundary  layers  on  electrode  surfaces  have  received  very  little  attention, 
reflecting  the  generally  unsatisfactory  state  of  our  understanding  of  electrical 
discharges  in  continuum  flows.  In  one  sense  this  lack  of  concern  on  the  part 
of  aerodynamicists  for  the  problem  of  current  emission  is  justified,  for  the 
plasma  sheath  which  forms  at  bounding  surfaces  is  not  a  boundary  layer  in 
the  ordinary  sense  of  fluid  dynamics,  but  is  more  aptly  described  as  a  problem 
in  surface  electro-chemistry.  Nevertheless,  it  is  to  be  hoped  that  this  problem 
will  not  continue  to  be  ignored. 

The  major  part  of  this  paper  will  be  devoted  to  a  discussion  of  MHD 
boundary  layers.  Section  11  deals  with  the  scaling  laws  appropriate  to  viscous 
MHD  flows,  while  Section  III  is  concerned  with  the  basic  mechanism  of 
diffusion  of  vorticity  and  current.  Section  IV  is  a  summary  of  present  under¬ 
standing  concerning  the  principal  types  of  MHD  boundary  layer  problems, 
including  turbulent  boundary  layers.  In  Section  V  we  discuss  some  plasma 
boundary  layers  in  which  magnetic  effects  are  absent. 


II.  Scaling  Laws 

The  interaction  of  a  conducting  fluid  with  an  electromagnetic  field 
must  be  characterized  by  a  larger  number  of  dimensionless  parameters 
than  is  needed  for  an  ordinary  fluid  flow.  The  bewildering  variety  of  these 
new  parameters  is  the  result,  not  of  an  attempt  to  inflate  the  currency  of 
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dimensional  analysis,  but  to  select  those  groups  which  most  aptly  describe 
the  principal  effects  in  a  particular  problem.  As  an  example,  consider  the 
conservation  of  mass  and  momentum  of  an  incompressible,  viscous,  electrically 
conducting  fluid: 


V 


V  =  0 


(2.  1) 


*1-}  Y--  ' 


p+jxB+pi'V  V 


(2.2) 


together  with  Maxwell's  equations  for  the  electromagnetic  field: 


V  X  E  = 


d  B 


(2.3) 


V  X  B  =  (2. 4) 

and  Ohm's  law  for  the  conduction  of  charge: 

j  =  CT  (E  +  y  X  B)  (2.  5) 

Following  the  usual  arguments,  the  viscous  Reynolds  number,  s  VL/y,  is 
the  ratio  of  inertia  force  to  viscous  force  in  a  steady  flow  as  determined  from 


Rationalized  mks  units  will  be  used  throughout. 
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the  second  and  last  terms  of  Eq.  (2.  Z).  By  a  similar  argument,  the  ratio 
of  magnetic  force  (j  x  B)  to  inertia  force  is  found  to  be  a  B^L/ p  V  s  S  when 
Eq.  (2.  5)  is  substituted  into  Eq.  (2.  2).  S  is  often  termed  the  magnetic 
interaction  parameter.  Alternatively,  one  might  compare  the  ratio  of 
magnetic  force  to  viscous  force,  thereby  defining  the  Hartmann  number, 

Mjj  =  LB  ^cf/p  V  .  Naturally,  R^,  S  and  are  not  independent,  but  are 
related  by  S. 

In  an  inviscid  flow  the  interaction  parameter,  S,  measures  the  effect 
of  the  magnetic  forces  on  the  flow.  For  an  MHD  boundary  layer,  however, 
it  would  seem  more  appropriate  to  compare  the  relative  effects  of  the  magnetic 
and  viscous  forces,  in  which  case  the  Hartmann  number  is  more  appropriate, 
especially  for  fully-devcloped  laminar  channel  flow  for  which  the  inertia  force 
is  zero.  Nevertheless,  in  some  boundary  layers  where  the  free  stream  flow 
disturbance  is  measured  by  S,  it  is  customary  to  retain  this  parameter  as 
descriptive  of  the  boundary  layer  flow  as  well. 

In  inviscid  flows  it  is  often  important  to  distinguish  between  the  component 
(B^)  of  magnetic  field  caused  by  currents  external  to  the  flow  and  the  component 
(B*)  induced  by  currents  flowing  within  the  fluid.  If  one  retains  this  distinction 
and  combines  Eqs.  (2.  4)  and  (2.  5)  to  write  Ohm's  law  as: 

V  X  B'  =  a  I^E  +  y  X  (B'  +  Bjj)J  (2.  6) 

then  the  condition  that  B'  be  much  less  than  B^  (and  thus  that  B^  should  appear 
in  the  definitions  of  S  and  is  that  cr  VL  =  magnetic 

Reynolds  number,  thereby  describes  the  effect  of  induced  currents  on 

the  applied  magnetic  field,  and  specifies  the  condition  under  which  the  magnetic 


field  may  be  considered  to  be  only  slightly  perturbed,  irrespective  of 
whether  or  not  the  flow  is  viscous. 

With  respect  to  MHlJ  boundary  layers,  it  can  be  expected  that  three 
independent  parameters,  such  as  R^,  and  S,  will  be  pertinent  in  specifying 

the  nature  of  the  flov».  In  some  degenerate  cases,  such  as  the  semi-infinite 
flut  plate  or  the  Rayleigh  problem  where  no  characteristic  length,  L,  exists, 
one  must  use  the  two  ratios: 


e 


(2.7) 


(2.8) 


€  ,  called  the  magnetic  Prandtl  number,  is  the  ratio  of  the  viscous  diffusivity 
[v)  to  the  magnetic  diffusivity  (K  ^  )”^)  ^nd  is  the  ratio  of  the  flow 

velocity  to  the  Alfven  velocity  ).  It  will  be  seen  below  in 

Section  III  that  both  of  these  parameters  arise  naturally  when  the  MHD  equations 
are  recast  in  the  form  of  diffusion  equations,  and  that  they  can  be  interpreted 
physically  as  the  ratio  of  diffusion  times  and  flow  times,  respectively. 

In  all  the  preceding  we  have  ignored  the  effect  of  compressibility, 
whose  existence  would  introduce  another  degree  of  freedom  and  the  usual 
parameter,  Mach  number.  With  the  benefit  of  hindsight,  it  appears  that  the 
effect  of  compressibility  within  the  boundary  layer  is  not  greatly  different  from 
what  is  ordinarily  the  case,  and  does  not  greatly  alter  our  previous  conclusions 
concerning  the  effect  of  the  magnetic  field. 
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We  have  patently  ignored  the  electric  field  in  our  order-of-magnitucie 
analysis,  and  it  plays  no  part  in  the  dimensionless  parameters  so  far 
considered  ^*.i8  is  a  result  of  the  fact  that  E  may  be  eliminated  between 
Eqs.  (2.  3]  and  (2.  5)  and  will  only  influence  the  flow  through  conditions  at 
the  boundary.  For  example,  in  channel  flows  E  appears  in  the  "loading 
parameter",  E/VB,  and  in  other  problems  determines  the  boundary  values  of 

aB 

1°^  -rr-  • 

The  point  of  view  so  far  maintained  is  one  which  describes  the  flow  in 
terms  of  continuum  conservation  laws  and  phenomenological  transport 
effects.  It  is  equally  important  to  delineate  the  limitations  imposed  by  the 
microscopic  structure  of  the  plasma.  To  do  this  without  an  extensive 
discussion  of  kinetic  tneory  is  difficult,  and  we  shall  be  content  to  emphasize 
only  the  salient  points  which  arc  most  pertinent  to  boundary  layers. 

If  we  adopt  as  the  basic  microscopic  length  the  Coulomb  radius,  r,: 


r  s 


2 


(2.  9) 


that  is,  the  distance  of  separation  of  two  atomic  charges  of  potential  energy 
equal  to  kT,  then  the  mean  free  path,^^,  of  the  charged  particles  is: 


^  =  (nrV' 


(2.  10) 
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except  for  a  logarithmic  factor  not  ordinarily  much  different  from  unity. 

For  a  continuum  description  of  the  plasma  motion,  it  is  certainly  necessary 
that  L  . 

The  lack  of  distinction  between  the  conduction  current  in  the  Ohm's 
law  of  Eq.  (Z.  5)  and  the  total  current  f  Eq.  (2.  4)  implies  the  usual 
assumption  that  the  charge  density,  p^,  is  so  small  that  the  convection 
current  can  be  neglected,  tha*-  is, 

(2.11) 

For  this  to  be  so,  the  geometric  scale,  L»  should  be  greater  than  the  Debye 
distance,  h,  which  is  the  distance  over  which  the  maximum  possible  charge 
density  (n^  e)  may  cause  a  potential  difference  kT/e: 


=^4  TT  r  n  )  ^^2 


By  combination  of  Eqs.  (2.  12)  and  (2.  10)  v/e  find  that 
3 


h 


n  r 
e 

4  t: 


{2.  12) 


1/2 


(2.  13) 


that  is,  the  Debye  distance  is  always  less  than  a  mean  free  path  whenever 
the  plasma  can  be  considered  to  be  a  perfect  gas  (volume  per  molecule  much 

3 

greater  than  x  ).  Thus,  charge  separation  imposes  a  less  stringent  condition 
than  that  of  the  continuum  assumption. 

At  electrode  surfaces  there  may  be  considerably  more  energy  than  kT 
available  to  a  charged  particle.  The  sheath  which  forms  at  such  surfaces  will 
have  a  thickness  comparable  to  that  obtained  by  replacing  kT/e  in  Eq.  (2.  12) 
by  the  electrode  potential  drop,  provided  this  does  not  exceed  the  mean  free 
path.  Except  in  non- continuum  flows,  these  sheaths  will  not  affect  the  fluid 
motion  within  the  boundary  layer. 

It  is  possible  to  replace  the  transport  coefficients  u  and  o’  in  the 
expressions  for  R,  S  and  Mj^  by  their  values  from  mean  free  path  theory 


CT  “ 


(2.  14) 


so  that  these  parameters  contain  explicitly  a  microscopic  length  in  a  fashion 
which  emphasizes  the  microscopic  structure  of  a  plasma.  More  to  the  point, 
however,  is  the  modification  to  Ohm's  law  (Eq.  2.  5)  which  results  from  the 
fact  that  the  electron,  ion  and  neutral  particle  mean  velocities  (  V  ,  V.  and  V  , 

'  -,e  -i»i  ~a 

respectively)  may  be  appreciably  different.  Momentum  gain  by  the  electrons 
due  to  their  electric  field  may  be  equated  to  their  momentum  loss  due  to 
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collisions  with  the  ions: 


-  n  e  (E  +  V  X  B)  ’ 
e  '  ^  ~e 


"e  "  Ji) 


Now  the  body  force,  j  x  B,  must  be  applied  to  the  heavy  particles,  and  a 

fraction,  n  /(n.  +  n  ),  of  this  force  must  be  transmitted  from  the  ions  to  the 
’  a' '  i  a* 

neutral  atoms  by  collisions  between  the  latter: 


(j  X  B)  = 


n.  m.  (V.  ~  V  } 
1  :  '-1  ^  a' 


If  the  mean  velocity  of  the  heavy  particles  is  defined  by: 


then  the  revised  Ohm’s  law  obtained  by  combining  Eqs.  (Z.  15),  (Z.  16)  and 
(2.  17)  is: 


u  r 

j=  o-(E  +  VxB)-  --g--  jxB 


<  “e  ■"e’  '  “i  ’’ia) 


(j  X  B)  X  B 
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The  second  term  on  the  right  side  is  the  Hall  field  which  results  from  the 
finite  value  of  the  electron-ion  collision  frequency  compared  with  the  electron 
cyclotron  frequency,  The  third  term  is  an  ion  slip  field  caused  by  the 

relative  motion  of  the  ions  and  atoms,  but  is  absent  when  the  plasma  is 
completely  ionized  (n^  =  0),  For  a  slightly  ionized  gas,  the  ion  slip  field 
is  comparable  to  the  Hall  field  when  o).  is  approximately  one. 

In  our  discussion  so  far  we  have  ignored  energy  conservation.  For  an 
incompressible  fluid,  the  variation  of  temperature  in  the  boundary  layer  will 
affect  the  flow  pattern  only  through  changes  in  the  transport  coefficients.  The 
coupling  in  a  compressible  flow  is  more  direct  because  of  the  additional  varia¬ 
tion  in  density.  Nevertheless,  the  principal  new  effect  on  the  macroscopic 
scale  would  by  the  joule  heating  (j  /a  per  unit  volume)  which  must  be  conducted 

to  the  wall  or  the  free  stream.  One  may  compare  the  ratio  of  the  joule  heating 
2 

to  the  heat  diffusion  (XV  T)  terms  in  the  energy  equation,  obtaining 


joule  heating 
heat  diffusion 


{Mj^  •  M)^ 


(2.  19) 


in  which  the  Prandtl  number  (Pj.)  and  Mach  number  (M)  are  introduced  as  is 
to  be  expected  in  compressible,  heat-conducting  flows. 

It  is  not  our  purpose  to  discuss  the  transport  properties  of  plasmas, 
except  insofar  as  they  will  affect  the  behavior  of  boundary  layers.  It  is 
sufficient  to  point  out  that,  as  a  gas  becomes  ionized,  the  thermal  conductivity 
increases  due  to  the  presence  of  the  light  electrons  while  the  viscosity  decreases 
for  high  degrees  of  ionization  because  the  Coulomb  cross  section  is  much  larger 
than  the  neutral  atom  cross  section.  Thus,  the  thermal  boundary  layer  thickness 
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is  much  greater  than  the  viscous  layer  thickness  in  a  completely  ionized  gas. 

On  the  microscopic  scale  we  can  again  expect  the  redistribution  of 
energy  to  modify  the  classical  situation  in  a  manner  similar  to  the  modifica¬ 
tion  of  Ohm's  law  occasioned  by  the  redistribution  of  momentum  between  the 
ions,  electrons  and  neutral  particles.  In  particular,  the  joule  heating  is  an 
energy  added  to  the  electrons  which  must  be  shared  quickly  with  the  particles 
with  which  they  collide  if  the  electron  temperature  (T^)  and  heavy  particle 
temperature  (T^)  are  to  be  very  nearly  the  same.  If  is  the  mean  time 
for  electrons  to  lose  their  thermal  energy  then 


or 


A 


(  w  %) 


(2.  20) 


If  energy  is  lost  by  elastic  collisions  of  irequency  \  then  T£=m.r^/m^ 
beC'JLUse  of  the  disparity  in  mass,  and  non-equilibrium  temperatures  will 
exist  only  if  Hall  effects  are  also  important.  If  inelastic  collisions  are 
effective,  then  usually  <  m.  T^/m^  and  temperature  equilibrium  is  more 
likely  to  be  obtained.  Whether  or  not  a  two-temperature  model  must  be  used 
within  the  boundary  layer  therefore  depends  quite  critically  upon  the  magnetic 
field  strength  (through  r^)  and  the  ratio  of  elastic  to  inelastic  electron 
cross  sections.  Thus,  one  additional  microscopic  parameter  must  be 
introduced  under  these  circumstances. 

As  is  the  case  for  boundary  layers  in  dissociated  gases,  cooled  plasma 
boundary  layers  will  recombine  ions  and  electrons  either  within  the  boundary 
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layer  by  gas  phase  recombination,,  or  at  the  wall  by  surface  recombination. 
The  diffusion  of  ion-electron  pairs  towards  the  wall  proceeds  by  ambipolar 
diffusion  in  which  the  high  diffusivity  of  the  electrons  induces  an  electric 
field  which  doubles  the  inherent  diffusivity  of  the  ions  through  the  neutral 
gas  in  the  vicinity  of  the  wall.  Since  the  electrical  conductivity  of  the 
boundary  layer  gas  is  very  sensitive  to  the  degree  of  ionization  when  the 
latter  is  less  than  a  few  percent,  there  may  again  be  strong  coupling  between 
the  effects  of  diffusion  and  recombination  on  one  hand  and  the  distribution  of 
electric  current  within  the  boundary  layer  on  the  other.  It  is  possible  that 
the  joule  heating  can  produce  ionization  (and  hence  a  source  of  diffusing 
ion-electron  pairs)  by  purely  equilibrium  thermal  effects  or  by  virtue  of 
superheating  of  the  electrons.  (We  will  leave  it  to  the  reader  to  define  the 
appropriate  parameter  for  the  ratio  of  ionization  source  to  ambipolar 
diffusion).  On  the  whole,  one  may  expect  that  many  of  the  techniques  used 
in  the  dissociating  gas  boundary  layer  may  be  adopted  for  this  case  as  well. 

In  the  foregoing  we  have  attempted  to  show  how  MHD  boundary  layers 
involve  several  additional  similarity  parameters  S,  .  .  .  )  which  can 

be  incorporated  in  the  usual  manner  in  the  continuum  analysis,  and  also 
microscopic  parameters  (ojt,  « •  •  )  which  usually  necessitate  a 

reformulation  of  the  classical  statements  such  as  Ohm’s  law,  energy  con¬ 
servation,  etc.  Within  this  framework  it  is  possible  to  treat  the  extremes 
of  liquid  metal  flow  in  a  tube  and  flow  of  a  partially  ionized  gas  through  a 
strong  magnetic  field  adjacent  to  a  cold  wall  which  "de-ionizes"  the  gas 
by  wall  recombination.  While  the  former  is  more  amenable  to  analysis, 
the  latter  is  of  great  practical  importance  and  must  undoubtedly  be  under¬ 
stood  more  thoroughly  than  is  presently  the  case.  Most  of  the  work 
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reviewed  below  is  much  more  closely  related  to  the  former,  rather  than 

the  latter,  problem. 

III.  The  Diffusion  of  Vorticity  and  Current 

One  way  to  characterize  the  difference  between  an  ordinary  boundary 
layer  and  an  MHD  boundary  layer  is  to  compare  the  processes  of  diffusion 
of  momentum  and  magnetic  field  energy.  To  simplify  this  analysis,  let 
us  consider  a  two-dimensional  incompressible  flow  with  co-planar  velocity 
and  magnetic  fields.  The  current,  vorticity,  and  electric  field  (neglecting 
H2dl  effects)  will  be  normal  to  the  plane  containing  V  and  B.  By  taking  the 
curl  of  Eqs.  (Z.  Z),  (Z.  4)  and  (2,  5),  we  find: 


3 

hr- 


+  V  ■  7  -uV 


(B  -  V)  j 


(3.1) 


V-  V  -  KV^I  j  =  -i-  (B-V)u  +-|-  V®  (B-V)x  V  (3.3)* 

The  lefthand  side  of  either  equation  expresses  the  usual  balance  of  diffusion 
and  convection  along  streamlines,  while  the  righthand  side  is  a  source  term 
which  is  proportionate  to  the  convection  of  the  complementary  quantity 
(j  or  w  )  in  the  direction  of  B,  It  is  this  convection  of  vorticity  or  current 
along  the  magnetic  field  lines  which  distinguishes  the  MHD  from  an  ordinary  flow, 

jjt  ^  1  Y 

The  last  term,  which  also  equals  -  Z/z^  V  (V*  V  )  x  B  ,  is  identically  zero 
for  the  Alfven  wave  and  Hartmann  boundary  layer  discussed  below,  and  also 
if  V  is  proportional  to  B.  For  small  Rj,^  and  constant  B^  it  may  be  neglected. 

In  the  following  discussion  it  will  be  considered  negligible. 
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A  possible  solution  to  Eqs.  (3.  1)  and  (3.  2)  for  tha  special  case 
of  =  K  is  one  for  which J  ^  *^^7  then 

be  written  as: 


d 

Tt 


+  (V  +  Ya^ 


<7 


(3.  3) 


which  expresses  the  propagation  of  a  diffusing  disturbance  with  a  velocity 
V  +  and  is  therefore  an  Alfven  wave.  In  MHD  boundary  layers  the 

boundary  conditions  on  j  and  w  generally  do  not  permit  their  being  every¬ 
where  proportional  in  the  ratio  p/p^  *  so  that  there  will  be  more  disturbances 
than  Alfven  waves  alone.  Nevertheless,  we  mast  expect  that  Alfven  waves  can 
be  excited  within  the  boundary  layer  and  thereby  carry  away  energy  and 
momentum. 

When  K  /  V  t  then  ^ j  /  u>  everywhere  within  the  Alfven 

wave,  but  the  current  and  vorticity  are  integrally  related  by: 


(3.  4) 


which  is  equivalent  to  the  statement  that  the  ratio  of  the  velocity  change 
across  the  wave  to  the  magnetic  field  change  (both  changes  being  parallel 
to  the  front)  equals  V^/B.  Thus,  by  Eq.  (3.  4)  the  ratio  of  total  vorticity 
to  total  current  in  an  Alfven  wave  is  • 

From  an  aerodynamic  point  of  view,  the  propagation  of  vortex  and 
current  elements  along  magnetic  field  lines  can  be  understood  by  noting 
that  the  force  —  V  x  w  on  a  vortex  convecting  with  a  velocity  V  can  be 


-14- 


balanced  by  a  force  j  x  B  provided  V  =  and  j  =5  Vp  w  , 

that  is,  the  current  and  vorticity  are  coincident  and  move  with  the 

Alfven  velocity. 

An  interesting  example  of  the  diffusion  of  vorticity  and  current 
within  a  boundary  layer  is  that  of  the  Hartmann  boundary  layer  which 


forms  on  a  fiat  plate  having  B  perpendicular  to  the  plane  of  the  plate. 
Assuming  the  flow  to  be  steady  and  to  depend  only  upon  the  distance 
y  normal  to  the  plate,  Eqs.  (3.  1)  and  (3.  2)  become: 


dj 


(3.5) 


(3.  6) 


Because  of  the  symmetry  of  the  plane  problem, 


constant  and  u>  and  j 


are  given  by: 


W  a  exp 

1-  - 

y 

I 

/  p 

V 

u. 

II 

<, 

°^i 

K 

<u 

where  5^^ 


(3.8) 
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Thus,  the  Hartma^m  boundary  layer  thickness,  6^,  is  the  usual  diffusion 
layer  thickness  for  a  stationary  source  in  a  fluid  of  diffusivity  -^kV  flowing 
with  a  velocity  V^.  In  this  sense  the  vorticity  and  current  are  convected 
toward  the  wall  with  the  Alfven  velocity  and  diffuse  away  from  the  wall  with 
the  geometric  mean  diffusivity  . 

One  may  use  this  picture  to  understand  the  development  of  an  MHD 
boundary  layer  on  an  infinite  flat  plate  suddenly  set  into  motion  (Rayleigh 
problem).  In  the  non-magnetic  case,  the  vorticity  diffuses  outv/ards  forming 
a  boundary  layer  of  thickness  =  -y^  .  In  the  MHD  case,  there  is  a 
tendency  for  vorticity  to  convect  towards  and  away  from  the  wall  with  the 
Alfven  velocity  as  well  as  to  diffuse  away  from  the  wall.  After  some  time, 
a  portion  of  the  vorticity  originated  by  the  plate  motion  convects  awa>  from 
the  wall  in  an  Alfven  wave  while  another  portion  establishes  a  steady  state 
balance  between  convection  toward  the  wall  and  diffusion  away  from  the  wall, 
that  is,  form^  a  Hartmann  boundary  layer. 

Ihe  relative  strength  of  these  two  effects  may  be  estimated  by 
applying  the  rule  that  the  total  current  which  flows  in  the  Hartmann  boundary 
layer  must  return  in  the  Alfven  wave  (assuming  an  insulating  plate)- 


(3.  9) 


The  velocity  change  across  the  Alfven  wave  (A  U^)  is  related  to  the  velocity 
change  across  the  Hartmann  layer  (A  U^j)  by: 


A 

“S~n7 


%  «A 


(3.  10) 
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Using  Eqs.  (3.  4),  (3.  7)  and  {3.  9),  this  may  be  evaluated  as: 

r- 

^  =vr  (3.U, 

Thus,  for  small  magnetic  Prandtl  number  (e),  the  Alfven  wave  propagation 
may  be  neglected  when  studying  boundary  layers  on  surfaces  with  a  magnetic 
field  normal  to  the  surface. 

Eqs.  {3,  1)  and  (3.  2)  are  the  basic  boundary  layer  equations  for  plane 
flows.  They  are  sometimes  linearized  with  respect  to  the  magnetic  field 
by  assuming  B  •  V  to  be  •  V  when  «1,  in  which  case  they  reduce  to: 

+  +  Y  ■  ('"a  |r| 

For  K  »  V  (i.  e.  ,  e  <.<  1),  within  the  viscous  boundary  layer  Eq.  (3.  12)  may 

9  2 

be  simplified  by  ignoring  the  terms  ^  V  •  V  compared  to  .  Either 

or  both  of  these  simplifications  are  commonly  Uwed  in  the  literature. 

The  question  of  MHD  wakes  which  propagate  along  the  magnetic  field 
lines,  even  in  the  upstream  direction  if  >  V  ,  has  received  considerable 

attention.  In  some  cases  the  flow  is  considered  inviscid  (1/=  0)  and  the  dis¬ 
turbance  originates  from  currents  induced  at  the  surface  of  a  body  which 
disturbs  the  flow.  This  might  be  classified  as  an  inviscid  magnetic  boundaxy 
layer  which,  in  the  most  general  case,  will  be  the  source  of  Alfven  waves 
spreading  into  the  flow.  It  is  also  possible  to  consider  the  complementary 
case,  a  viscous,  superconducting  boundary  layer  (K  =  0)  in  which  the  roles 
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of  j  and  w  are  interchanged  in  accordance  with  the  symmetry  of  Eqs. 
(3.  I)  and  (3.2). 

For  a  plasma,  the  value  of  c  is  found  to  be 


where  n^  is  the  electron  density  in  particles  per  cubic  centimeter. 

Because  of  the  steep  temperature  dependence  of  both  v  and  a  ,  e  «  1 

5  16 

for  low  temperatures  (T  <  10  at  =  10  )  but  e  »  1  at  high  temperatures. 

However,  the  conditions  under  which  €  »  1  ff>r  a  plasma  generally  imply 
such  strong  magnetic  fields  that  the  continuum  MHD  assumptions  no  longer 
hold. 


IV.  MHD  Boundary  Layers 

The  MHD  Rayleigh  problem,  a  particular  case  of  which  was  discussed 

47 

above,  was  first  treated  by  Rossow  for  the  case  of  £  =  1.  However,  the 
boundary  conditions  on  the  electric  field  used  by  Rossow  were  not  physically 
meaningful,  leaving  arisen  through  a  misapprehension  concerning  "moving" 

49 

and  "stationary"  magnetic  fields.  A  proper  solution  for  a  superconducting 
plate  has  been  obtained  by  Chang  and  Yen^,  and  a  further  discussion  is  given 
by  ■Williams. 
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The  corresponding  Stokes  problem  (oscillating  infinite  plate,  B 

40 

perpendicular  to  the  plate),  was  solved  by  Ong  and  Nicholls  for  e  «  1 

21 

and  «  1.  More  complete  solutions  were  obtained  by  Hide  and  Roberts 
2 

and  Axford  for  other  limiting  cases. 

The  incompressible  flow  over  a  semi- infinite  flat  plate  with  B 

14  6 

parallel  to  the  flow  was  shown  by  Greenspan  and  Carrier  '  to  involve 

an  upstream  wake  if  V  <  so  that  steady  solutions  exist  under  such 

1 9 

circumstances  only  for  a  finite  plate.  Asymptotic  solutions  for  large 
or  small  e  were  given  by  Glauert.  For  this  configuration,  the  disturbance 
to  the  magnetic  field  is  caused  by  the  displacement  effect  of  the  viscous 
boundary  layer. 

The  magnetic  Blasius  problem  with  B  perpendicular  to  the  wall  has 

Q 

been  discussed  by  Fay  ,  who  pointed  out  the  significance  of  the  Hartmann 

41  58 

layer,  and  related  incompressible  flews  have  also  been  considered.  * 

4  44  45 

More  attention  has  been  given  to  variable  a  hypersonic  boundary  layers  *  * 

31  32  34  4 

and  compressible  boundary  layers.  '  '  In  the  former  case.  Bush 

found  two  possible  solutions  under  some  circumstances,  probably  due  to 

the  very  strong  dependence  of  <t  upon  the  temperature,  which  peaks  in  the 

center  of  the  hypersonic  boundary  layer  when  the  wall  is  cooled.  Kerrebrock^^* 
38 

and  Napolitano  and  Pozzi  have  discussed  flat  plate  boundary  layers  with 
B  parallel  to  the  wall  but  perpendicular  to  the  flow.  Kerrebrock^"^  also 
considered  the  case  of  a  current  flowing  into  the  fluid  from  the  plate  (electrode). 
Because  of  its  great  practical  interest,  the  MHD  boundary  layer  at  the 

stagnation  point  of  a  blunt  body  having  B  normal  to  the  surface  has  been  very 
46 

thoroughly  studied.  Rossow  showed  that  the  heat  transfer  to  a  cylinder 
would  be  reduced  by  the  magnetic  field.  Meyer^^  pointed  out  that  the 
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principal  effect  of  the  magnetic  field  was  to  alter  the  inviscid  flow, 
thereby  reducing  the  stagnation  point  velocity  gradient  andi  hence,  heat 

transfer.  A  discussion  of  the  inviscid  flow  calculation  and  summary  of 

Z4  ZS 

the  heat  transfer  effects  have  been  given  by  Kemp  *  ,  but  the  inviscid 

1 8 

flow  calculations  have  been  subject  to  some  debate.  Recent  developments 

23  39 

include  an  incompressible  solution,  extension  to  rarefied  gas  flow, 

57  5 

and  an  attempt  to  find  solutions  for  large  Bush  has  shown  that 

the  reduction  of  o*  within  a  cooled  stagnation  point  boundary  layer  causes 
an  overshoot  in  the  velocity  profile  and  a  lesser  reduction  in  heat  transfer 
than  would  be  the  case  for  constant  a.  This  overshoot,  which  is  caused 

by  the  fact  th<  t  the  pressure  gradient  at  the  stagnation  point  is  unaffected 

35 

by  the  magnetic  field,  has  been  discussed  by  Lykoudis. 

The  simolest  channel  flow  is  Couette  flow,  which  was  treated  by 

3  54 

Bleviss  for  a  variable  property  hypersonic  flow.  Tao  considered  the 

incompressible  time-dependent  Cottette  flow  which  is  related  to  the 

Rayleigh  problem.  Hartmann  type  channel  flows  have  been  investigated  by 

Shercliff^^  for  the  case  of  circular  pipes  and  by  Globe for  the  case  of  an 

annulus.  The  time-dependent  channel  flow  has  been  investigated  by  Bundgren, 

Atabek  and  Chang.  In  Che  3t*!ady  flow,  the  Hartmann  number  is  the 

primary  parameter  defining  the  flow  while  the  Reynolds  number  is  also 

involved  in  the  unsteady  flow. 

9 

Very  little  attention  has  been  given  to  Hall  effects.  Fay  has  shown 
that  the  Hall  current  causes  a  cross-flow  in  the  Hartmann  boundary  layer. 
Sutton  cind  Sherman^^  considered  Hall  effects  in  Poiseuille  flow  and  also 
found  the  s^e  cross-flow  due  to  Hall  currents.  The  principal  effect  of 
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permitting  Hall  currents  to  flow  appears  to  be  that  of  reducing  the  elec- 

2 

trical  conductivity  approximately  by  the  factor,  I  +  (w^  t^)  . 

Studies  of  MHD  boundary  layer  stability  have  shown,  in  general, 

that  the  magnetic  field  has  a  stabilizing  effect.  For  a  flat  plate  boundary 

43  1 

layer  with  B  parallel  to  the  flow,  Rossow  and  Arkhipov  found  that  the 
minimum  critical  Reynolds  number  more  than  doubles  for  an  interaction 
parameter  somewhat  less  than  unity.  For  the  corresponding  free  boundary 
layers,  which  are  always  unstable  in  the  absence  of  a  magnetic  field,  the 

flow  may  be  stabilized  for  a  finite  wave  number  by  a  sufficiently  ^arge 
8  13 

magnetic  field  and  is  absolutely  stable  for  all  R.y  if  £  >  1  and  V  <  V^. 

For  Poiseuille  flow  with  B  parallel  to  the  flow,  Stuart^^  found  a  sixfold 
increase  in  the  minimum  critical  Reynolds  number  for  S  =  0.  10.  When 
B  is  perpendicular  to  the  wail,  the  minimum  critical  Reynolds  number 

29 

also  increases  with  increasing  approaching  a  value  of  50,000 

when  Mj_j  »  1.  {The  experimental  value  of  at  which  transition  occurs 

was  found  by  Murgatroyd  to  be  about  250.  Since  Ry/Mj^  is  the  value  of  the 
viscous  Reynolds  number  based  on  the  Hartmann  thickness,  it  would  appear 
that  in  the  experiments  the  principal  effect  of  the  magnetic  field  is  to  alter 
the  basic  flow  rether  than  to  stabilize  it. ) 

Most  of  the  approaches  to  turbulent  MHD  flows  have  been  speculative 

in  nature,  being  patterned  after  classical  analyses  of  turbulence. 

28 

Kovasznay  has  investigated  the  resistivity  of  a  turbulent  plasma  in  which 
the  current  paths  are  greatly  lengthened  by  the  random  V  x  B  field, 
and  the  turbulent  energy  is  provided  by  random  j  x  B  forces.  It  v/as 
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suggested  by  L.ykoudis  that  transition  to  turbulence  in  channel  flows 
will  occur  when  the  Reynolds  number  of  the  Hartmann  layer  exceeds 
that  of  the  laminar  sub-layer  (about  200)  in  a  turbulent  flat-plate  boundary- 
layer.  This  hypothesis  is  in  reasonable  agreement  with  Murgatroyd's 
experiments, 

15-17 

The  general  discussions  of  wakes  have  been  concerned  with 

the  propagation  of  Alfven  waves  originating  from  the  surface  of  a  finite 

body,  and  make  use  of  the  magnetic  equivalent  of  the  Oseen  approximation. 
48 

Sears  has  discussed  magnetic  boundary  layers,  but  no  solutions  oi  such 

problems  for  R^  not  small  have  been  obtained.  MHD  jet  flow^^  and  natural 
42 

convection  flow  have  recently  been  studied,  at  least  within  the  boundary 
layer  approximation. 


V,  Electronic  Boundary  Layers 

Boundary  layers  in  plasmas  with  no  magnetic  field  present  arc  much 
more  nearly  like  ordinary  boundary  layers.  The  principal  interest  in  such 
boundary  layers  lies  in  the  effects  of  the  presence  of  electrons  on  the  thermal 

9 

conductivity,  which  might  considerably  enhance  the  heat  transfer.  An 
experimental  study  of  the  increase  in  heat  trcuisfer  due  to  the  passage  of  a 
current  through  such  a  boundary  layer  has  been  made  by  Fay  and  Hogan, 
who  found  that  the  inci  ease  at  either  cathode  or  anode  was  proportional  to 
the  current.  Talbot^*  has  analyzed  the  plasma  sheath  inside  a  stagnation 
point  boundary  layer  on  a  Langmuir  probe,  but  no  experimental  data  is 


available. 


VI.  Conclusion 


In  surveying  the  progress  to  date,  one  is  struck  by  the  comparatit’-e 
richness  of  the  theoretical  analyses  of  plasma  boundary  layers  as  compared 
to  the  paucity  of  experimental  evidence.  As  noted  previously,  this  reflects 
the  relatjve  case  of  extrapolating  existing  boundary  layer  techniques  as 
opposed  to  developing  an  entirely  new  set  of  laboratory  experiments  and 
techniques.  W’liile  one  might  be  tempted  to  conclude  tliai  our  theoretical 
understanding  of  the  problem  is  far  ahead  of  our  laboratory  esrperience,  it 
IS  more  likely  than  not  that  we  have  been  analyzing  the  wrong  problems.  For 
the  variety  of  physical  phenomena  exhibited  by  plasmas  is  certainly  no  less 
than  that  to  be  seen  in  ordinary  fluids,  and  quite  possibly  is  even  more 
extensive  and  of  greater  pra.ctical  utility'.  It  seems  highly  unlikely  .oat  an 
adequate  understanding  of  fluid  plasma  behatdor  can  be  manufactured  from 
the  ingredients  of  classical  fluid  meciianics  and  electromagnetic  theory  aione» 

In  addition  to  an  increased  emphasis  on  experimental  research,  one 
would  hope  for  a  better  otiderslanding  of  the  microscopic  processes  which 
play  such  an  important  role  in  transport  properties  and  the  tendency  toward 
maintenance  of  thermodynamic  equilibrium,  Beciuse  of  the  difficulty  of 
creating  and  maintaining  a  plasma,  one  is  inexatably  forced  to  work  under 
marginal  circumstances  which  make  necessary  a  quantitative  knowledge  of 
those  processes  by  which  energy  is  lost  and  the  plasma  degraded. 

Finally,  one  must  be  prepared  for  more  constructive  concepts  of 
plasma  dynamics  if  such  seem  useful.  Only  in  a  very  limited  sense  is  a 
plasma  an  anisotropic,  electrically-conducting  fluid.  It  might,  for  example, 
be  more  convenient  to  consider  it  to  be  a  syslen:i  of  hydromagnetic  waves. 


for  which  an  appropriate  kinetic  theory  and  macroscopic  conservation 


laws  must  be  devised.  Given  such  alternative  s»  the  problem  of  boundary 
effects  would  undoubtedly  acquire  a  different  significance. 
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